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Abstract 

A Monte Carlo event generator has been developed assuming thermal production of hadrons. 
The system under consideration is sampled grand canonically in the Boltzmann approximation. 
A re-weighting scheme is then introduced to account for conservation of charges (baryon number, 
strangeness, electric charge) and energy and momentum, effectively allowing for extrapolation of grand 
canonical results to the microcanonical limit. This method has two strong advantages compared to 
analytical approaches and standard microcanonical Monte Carlo techniques, in that it is capable of 
handling resonance decays as well as (very) large system sizes. 
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I. INTRODUCTION 

The statistical hadronization model, first introduced by Fermi |l| and Hagedorn j^, has been 
remarkably successful in the description of experimentally measured average hadron production 

ri ri ri 

yields in heavy ion collisions ranging from SIS (3[, and AGS over SPS M to RHIC m 
energies. Over time this has led to the establishment of the 'chemical freeze-out line' which 
is now a vital part of our understanding of the phase diagram of strongly interacting matter. 
Model predictions for the upcoming LHC and future FAIR j^, 9] experiments largely follow 
these trends. 

Somewhere above this freeze-out line in the phase diagram we expect, in general, a phase 
transition from hadronic degrees of freedom to a phase of deconfined quarks and gluons, gen- 
erally termed the quark gluon plasma; and more specifically, a first order phase transition at 
low temperature and high baryon chemical potential, and a cross-over at high temperature and 



low baryon chemical potential. In 



emerge. For recent reviews see 
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between, a second order endpoint or a critical point might 
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Fluctuation and correlation observables are amongst the most promising candidates sug- 
gested to be suitable for signaling the : 
them. For recent reviews here see Il2l. 



brmation of new states of matter, and transitions between 



15fl. 

The statistical properties of a sample of events are, however, certainly not solely determined 
by critical phenomena. More broadly speaking, they depend strongly on the way events are 
chosen for the analysis, and on the information available about the system. 

The ideal gas approximation of the statistical hadronization model will again serve as our 
testbed. Its strong advantage is that it is simple, and to some extent intuitive. Given its success 
in describing experimentally measured average hadron yields, and its ability to reproduce low 
temperature lattice susceptibilities |l6|, the question arises as to whether fluctuation and cor- 
relation observables also follow its main line. Critical phenomena (and many more), however, 
remain beyond the present study. 

Conventionally in statistical mechanics three standard ensembles are discussed; the micro- 
canonical ensemble (MCE), the canonical ensemble (CE), and the grand canonical ensemble 
(GCE). In the MCE'^ one considers an ensemble of microstates with exactly fixed values of 



^ The term MCE is also often applied to ensembles with energy but not momentum conservation. 
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extensive conserved quantities (energy, momentum, electric charge, etc.), with 'a priori equal 



probabilities' of all microstates (see e.g. [17l]). The CE introduces the concept of temperature 
by introduction of an infinite thermal bath, which can exchange energy (and momentum) with 
the system. The GCE introduces further chemical potentials by attaching the system under 
consideration to an infinite charge bath^. Only if the experimentally accessible system is just a 
small fraction of the total, and all parts have had the opportunity to mutually equilibrate, can 
the appropriate ensemble be the grand canonical ensemble. 

A statistical hadronization model Monte Carlo event generator affords us with the possibility 
of studying fluctuation and correlation observables in equilibrium systems. Data analysis can 
be done in close relation to experimental analysis techniques. Imposing global constraints 
on a sample is always technically a bit more challenging. Direct sampling of MCE events (or 
microstates) has only been done in the non-relativistic limit 13]. Sample and reject procedures, 
suitable for relativistic systems, become rapidly inefficient with large system size. However, they 



have the advantage of being very successful for small system sizes |l9l . |20[ |. 

In this article we try a different approach: we sample the GCE, then re-weight events ac- 
cording to their values of extensive quantities, and approach the sample-reject limit (MCE) 
in a controlled manner. In this way one can study the statistical properties of a global equi- 
librium system in their dependence on the size of their thermodynamic bath. As any of the 
three standard ensembles remain idealizations of physical systems, one might find intermediate 
ensembles to be of phenomenological interest too. 

We study the first and, in particular, second moments of joint distributions of extensive quan- 
tities. We concentrate mainly on particle number distributions and distributions of 'conserved' 
charges, and discuss the infiuence of acceptance cuts in momentum space, conservations laws, 
and resonance decay on the statistical properties of a sample of hadron resonance gas model 



events. We extend our previous studies of ideal particle and anti-particle gases [2l|, l22| and of 
gases of altogether massless particles 



The numerical code has been written for inclusion into the already existing THERMUS 



package [2J]. We make frequent use of the functionality provided by the ROOT framework [25 1. 

The paper is organized as follows: In Section [III the basic ideas of this article are formu- 
lated. The GCE Monte Carlo sampling procedure is described in Section IIIII The first and 



Note that a system with many charges can have some charges described via the CE and others via the GCE. 
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second moments of the distributions of fully phase space integrated extensive quantities are 
then extrapolated to the microcanonical limit in Section IIVI Section |V] contains an analysis of 
GCE momentum spectra. The momentum space dependence of correlations between conserved 
charges is studied in Section IVII Section IVIII then deals with multiplicity fluctuations and 
correlations in limited acceptance and their extrapolation to the MCE limit. A summary is 
given in Section IVIIIi 

II. STATISTICAL ENSEMBLES WITH FINITE BATH 



We start out as Patriha [17|], and Challa and Hetherington (26|, but quickly take a different 
route. 



Let us define two microcanonical partition functions, i.e. the number of microstates, for 
two separate systems. The first system is assumed to be enclosed in a volume Vi and to have 
fixed values of extensive quantities Pf = {Ei, P^^i, Py^i, P^^i), and Q{ = {Bi, Si,Qi), while 
the second system is enclosed in a volume V2 and has fixed values of extensive quantities 
P2 = (£2, Px,2, Py,2, Pz,2), and Q2 = (-82,52,(52), where E is the energy of the system, Px,y,z 
are the components of its three-momentum, and B, S, and Q, are baryon number, strangeness 
and electric charge, respectively. Thus we have: 

Z(\/i,Pf,Qi) = 5^ Z^.(yi,Pf,Qi) , and Z {V2, P2^ , Qi) , (1) 

{Nl} 

where Zj^i{Vi,Pi,Q{) denotes the number of microstates of system 1 with additionally fixed 
multiplicities Nl of particles of species i. Suppose that system 1 and system 2 are subject to 
the following constraints: 

Vg = V, + V2, (2) 

Pg^ = P^ + P^ , (3) 

Qi = Qi + Qi ■ (4) 

We can then construct the partition function Z{Vg, P^", Q^g) of the joint system as the sums over 
all possible charge and energy-momentum split-ups: 

ZiVg, P^, Ql) = J2Y1 ^(^^ - ^9 - A", QI - QI) ZiVr, Pf, Q{) . (5) 
iPi} {QU 
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Next we construct the distribution of extensive quantities in the subsystem Vi. This is given by 
the ratio of the number of all microstates consistent with a given charge and energy-momentum 
split-up and a given set of particle multiplicities to the number of all possible configurations: 

We then define the weight factor W{Vi, P^, Q{; Vg, Pj^, Q^) such that: 

p(pf,gi,ivi) = w{VuP,',QhVg,Pg^,Qi) z^.(Vi,pf,gi) . (7) 

By construction, the first moment of the weight factor is equal to unity: 

w = EEE w{VuPr,Qi;Vg,p^^,Qi) z^.{v^,pr,Q{) = i, (s) 

{pn{Q{}m} 

as the distribution is properly normalized. 

The weight factor , Q{] Vg, P^, generates an ensemble with statistical properties 

different from the limiting cases Vg Vi (MCE), and Vg ^ oo (GCE). This effectively allows 
for extrapolation of GCE results to the MCE limit. In the thermodynamic limit (Vi sufficiently 
large) a family of thermodynamically equivalent (same densities) ensembles is generated. In 
principle any other (arbitrary) choice of Vr(Vi, Pf , Q{] Vg, P^l", Q^) could be taken. In this work 
we confine ourselves, however, to the situation discussed above. Please note that all microstates 
consistent with the same set of extensive quantities {Pi, Q{) have 'a priori equal probabilities'. 

In the large volume limit, ensembles are equivalent in the sense that densities are the same. 
The ensembles defined by Eq. ([7]) and later on by Eq. ffTTl) are no exception. If both Vi and Vg 
are sufficiently large, then the average densities in both systems will be the same, Q-'g/Vg and 
Pg /Vg respectively. The system in Vi will hence carry on average a certain fraction: 

A ^ Vl/Vg , (9) 

of the total charge and four-momentum P^, i.e.: 

{Q{) = XQl, and (Pf) = X Pg^ . (10) 

By varying the ratio A = Vi/Vg, while keeping {Q{) and (P^) constant, we can thus study 
a class of systems with the same average charge content and four-momentum, but different 
statistical properties. 
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A. Introducing the Monte Carlo Weight W 

Since Eq.([7]) poses a formidable challenge, both mathematically and numerically, we write 
instead: 

where the distribution of extensive quantities Pi, Q\ and particle multiplicities A'^^ of a GCE 
system with temperature T = volume Vi, chemical potentials /ij and collective four- velocity 
is given by: 

p,,e(Pr,Ql,ivi|/?,n„/i,) ^ z{yl,f^,u,,^^,) ^Niiyi,PlQ{) , (12) 

where /i-,- = {fiB, fJ'S, fJ'Q), summarizes the chemical potentials associated with baryon number, 
strangeness and electric charge in a vector. The normalization in Eq. (fT2|) is given by the 
GCE partition function u^, /ij), i.e. the number of all microstates averaged over the 

Boltzmann weights e^-^i'"''^ and e'^i'^-'^: 

Z{Vi,P,u„f^,) = E E E ^"^'"'^ Z^.{Vi,Pr,Q\) ■ (13) 

{Pi} {Oil {^J} 

The new weight factor W^i'''^i'^3''*^s(Vi; u^, /ij) now reads: 

„ , „ , p-iP^-PtW(^ piQi-Qi)HP 

^ Z{V,-Vi,P^-P^,Ql-Q{) ^4) 

z{v,^pi:.Qi) 

In the case of an ideal (non- interacting) gas, Eq. (fT4l) can be written {21! . I27I as: 

Z^9.«s(\/^,/5,u^,/ij) 

The advantage of Eq. (fTT|) . compared to Eq.Q, is that the distribution Pgce{Pi, Qi, Nl\P, u^, fij) 
can easily be sampled for Boltzmann particles, while a suitable approximation for the weight 
}^Pi^,Qi;P,^,Qi(^Vi; Vg\(3, u^, jjj) is available. 

Again, by construction, the first moment of the new weight factor is equal to unity: 

W = E E E p^^^(pM,gj,iVl|/5,u„/i,) = 1 . (16) 

{Pi^}{Qi}m} 



In principle, Eq.(I7]) and Eq.l jTTl) are equivalent. In fact, Eq.Q can be obtained by taking 
the limit {hb, jJ^Si I^q) = (0,0,0), = (1,0,0,0), and /3 — > of Eq.(|TT]). However, as one can 
already see, (W^) 7^ {W^). Higher, and in particular the second, moments of the weight factors 
W and W are a measure of the statistical error to be expected for a finite sample of events. 
The larger the higher moments of the weight factor, the larger the statistical error, and the 
slower the convergence with sample size. Please see also Appendices \M and El 

As GCE and MCE densities are the same in the system V^, these values are effectively 
regulated by intensive parameters /5, /ij and n^. In essence, if you want to study a system with 
average (Qi), then sample the GCE with (Qi) and calculate the weight according to Eq. (fT5|l . 
This will result in a low statistical error for finite samples (as shown in later sections), and 
allow for extrapolation to the MCE limit. 

We will now first calculate the weight factor Eq.f llSI) and then take the appropriate limits. 
With the appropriate choice of /5, /Xj and the calculation of Eq. (fT5i) is particularly easy in 
the large volume limit 27l |. 



B. Calculating the Monte Carlo Weight W 

In this article, the total number of (potentially) conserved extensive quantities in a hadron 
resonance gas isL = J + 4 = 3 + 4 = 7, where J = 3 is the number of charges (5, S, Q) 
and there are four components of the four-momentum. Including all extensive quantities into 
a single vector: 

= {Q\ Pn = (B, S, Q, E, P., Py, P,) , (17) 
the weight Eq. (fT5|) can be expressed as: 

The general expression for the partition function Z^' {V, P,u^, fXj) in the large volume limit 



reads 



(19) 



where: 

e = {Q'-Vk'i) {a-')[ , (20) 
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and: 



1/2 



Here ki and K2 are the GCE vector of mean values and the GCE covariance matrix respectively. 
The values of j3, fij and are chosen such that: 



dQ} 



0,. 
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The approximation (fT9l) gives then a reliable description of Z^'^a around the equilibrium value 
2g = ^g^i^ provided Vg is sufficiently large. The charge vector, Eq. (l20|) . is then equal to the 



null- vector = 0; (Q^ = VgK[). 

For the normalization in Eq. (fT8|) we then find: 



Ql =Ql 



Z{Vg,P,U^,fij) 



exp(O) . 



(23) 



For the numerator we obtain: 



Q' =Q' 



Z{Vg -Vl,(3,Uf„fIj) ( \ 1 

exp 



{^Txiyg - Vx)YI^ deter 
where in Eq. (1241) we write for the charge vector Eq. (l20l) : 



2 {Vg - 



i% \ , (24) 



Then, using Qg = Qg^^q = ^9^1' 

{AQ2)' = iQg-Qi)'-{Vg-Vi)^1 = -{Q,-Vim)' 
Substituting Eq.([23D and Eq.([24D into Eq.([T8]) yields: 



(25) 



(26) 



W^-^a(^V^■,Vg\(3,u„^l,) 



0' =o' 



{2nVgfl'^ deter 



(2vr {Vg - V,)Y/^ det a ""^H 2 {Vg - V,) ^^^^ 
The GCE partition functions are multiplicative in the sense that Z{Vi, l3,u^, fij) Z{Vg — 
Vi, P,Ufj,, fij) = Z{Vg, P,u^, ^j), and thus the first term in Eq. (!27|) is equal to unity. Now 
using Eq.(l9]), A = Vi/Vg, we can re- write Eq.( l27I) as: 



wQ-Qa{Vi;Vg\p,u,,f,,: 



o' =0' 
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Model parameters are hence the intensive variables inverse temperature /5, four-velocity 
and chemical potentials fi^ , which regulate energy and charge densities, and collective motion. 
Provided Vi is sufficiently large, we have defined a family of thermodynamically equivalent 
ensembles, which can now be studied in their dependence of fiuctuation and correlation ob- 
servables on the size of the bath V2 = Vg — Vi. Hence, we can test the sensitivity of such 
observables, for example, to globally applied conservation laws. The expectation values (...) 
are then identical to GCE expectation values, while higher moments will depend crucially on 
the choice of A. 

C. The Limits of W 

The largest weight is given to states for which = 0, i.e. with extensive quantities 
Qi = Qi eq ■ Hence, the maximal weight a microstate (or event) at a given value of A = Vi/Vg 
can assume is Wmax "(Vi; Vg\(3, m^, fij) = (1 — X)"^'"^. Taking the limits of Eq.( l28i) . it is easy to 
see that: 

hm W^'^■'^'^{V^■Vg\(3,u^,^J) = 1. (29) 

A— ^0 

I.e. for A = we sample the GCE, and all events have a weight equal to unity. Hence, we also 
find (W^) = 1 and therefore ((AW)^) = 0, implying a low statistical error. For A ^ 1, we 
effectively approach a "sample- reject" procedure, as (for instance) used in jl^l, and: 

lim W^'^'^HVi;Vg\P,u^,fi,) (X 5{Q[-Vik[). (30) 

However, as now not all events have equal weight, ((AW)^) grows and so too the statistical 
error of finite samples. Also, the larger the number L of extensive quantities considered for 
re-weighting, the larger will be the statistical uncertainty. 

III. THE GCE SAMPLING PROCEDURE 

The Monte Carlo sampling procedure for a GCE system in the Boltzmann approximation 
is now explained. The system to be sampled is assumed to be in an equilibrium state 
enclosed in a volume Vi with temperature T = j3~^ and chemical potentials fij = (yU^, /i^, /xq). 
Additionally, the system is assumed to be at rest. The four-velocity is then = (1, 0, 0, 0) and 
the four-temperature is = 0, 0, 0). In this case, multiplicity distributions are Poissonian, 
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while momentum spectra are of Boltzmann type. 



The GCE samphng process is composed of four steps, each discussed below. 



1. Multiplicity Generation 

In the first step, we randomly sample multiplicities of all particle species i considered 
in the model. The expectation value of the multiplicity of thermal Boltzmann particles in the 
GCE is given by: 

{N\) = 1^ T K, (^) e-/^ . (31) 

Multiplicities {N\}n are randomly generated for each event n according to Poissonians with 
mean values {Nl): 

P{Nl) = e"<^^) . (32) 

In the above, and Qi are the mass and degeneracy factor of a particle of species i respectively. 
The chemical potential fii = ^ijql = fisbi + fisSi + fiQQi, where = {bi,Si,qi) represents the 
quantum number content of a particle of species i. 



2. Momentum Spectra 

In the second step, we generate momenta for each particle according to a Boltzmann spec- 
trum. For a static thermal source spherical coordinates are convenient: 

M^'-^^' e-«/- . ,33) 
These momenta are then isotropically distributed in momentum space. Hence: 

Vx = bl sin 6 cos (p , (34) 

Vy = \p\ sin 9 sm(f) , (35) 

= \p\ cos 9 , (36) 

e = + , (37) 

where px, Py, and pz are the components of the three-momentum, e is the energy, and |j9| = 
■\/pI + Py+ vI is the total momentum. The polar and azimuthal angles are sampled according 
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to: 

9 = cos"^ [2 (x - 0.5)] , (38) 
(j) = 2 TT (x - 0.5) , (39) 

where x is uniformly distributed between and 1. Additionally, we calculate the transverse 
momentum px and rapidity y for each particle: 



pt = \ pI+pI , (40) 



^ = (41) 



2 Ve-P. 

Finally, we distribute particles homogeneously in a sphere of radius ri and calculate decay times 
based on the Breit-Wigner width of the resonances. 

3. Resonance Decay 

The third step (if applicable) is resonance decay. We follow the prescription used by the 



authors of the THERMINATOR package 29|, and perform only 2 and 3 body decays, while 
allowing for successive decay of unstable daughter particles. Only strong decays are considered, 
while weak and electromagnetic decays are omitted. Particle decay is first calculated in the 
parent's rest frame, with daughter momenta then boosted into the lab frame. Finally, decay 
positions are generated based on the parent's production point, momentum and life time. 
Throughout this article, always only the lightest states of the following baryons: 

p n A S+ S" S" S° n- (42) 

and mesons: 

7r+ TT- 7r° K- (43) 

are considered as stable. The system could now be given collective velocity u^. 

4- Re-weighting 

In the fourth step, we calculate the values of extensive quantities for the events gener- 
ated by iterating over the particle list of each event. For the values of extensive quantities 
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Q{ ,^ = {Bi^n, Si^n, Qi,n, -Ei,n, Px,i,n, Py,i,n, Pz,i,n) m subsystem Vi of eveiit n we write: 

particles i„ 

where q^^ = {bi„, Si^,qi,^,ei,^,p^^i„,Py^i„,P;,^iJ is the 'charge vector' of particle i in event n. Based 
on we calculate the weight w;„ for the event: 

= W^'---'^HVi;Vg\^,u^,fi,) , (45) 

according to Eq. (!28|) . Please note that all microstates with the same set of extensive quantities 
Q[ „ are still counted equally. 

IV. EXTRAPOLATING FULLY PHASE SPACE INTEGRATED QUANTITIES TO 
THE MCE 

We now attempt to extrapolate fully phase space integrated grand canonical results to 
the microcanonical limit. For this we iteratively generate, re- weight, and analyze samples of 
events for various values of \ = Vi/Vg. By construction of the weight factor W, Eq. (l28l) . 
we extrapolate in a systematic fashion such that, for instance, particle momentum spectra 
as well as mean values of extensive quantities remain unchanged. On the other hand, all 
variances and covariances of extensive quantities subject to re-weighting converge linearly to 
their microcanonical values. 

This can be seen from the form of the analytical approximation to the grand canonical 
distribution of (fully phase space integrated) extensive quantities Pgce(Qi) (from Eq. (fT9l) ): 

^-(Si) - p^^y^ -p (4 i rt,) . (46) 

where the variable is given by Eq. fl20|) . Now taking the weight factor Wa, Eq. fl28l) . [a and 
are the same in both equations) we obtain for the distribution P\{Q[) of extensive quantities 
Q[ in subsystem 1: 

Px{Q[) ^ P,UQ[) (47) 

This is essentially the same multivariate normal distribution as the grand canonical version 
Pgce{Qi), however hnearly contracted. We will compare Monte Carlo results to Eq. fj48|l . 
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The Monte Carlo output is essentially a distribution Pmc{Xi, X2, X^, ...) of a set of observ- 
ables Xi, X2, X3, etc. For all practical purposes this distribution is obtained by histograming 
all events n according to their values of Xi „, X2,„, X^^n, etc. and their weight Wn- One can 
then define moments of two observables Xj and Xj through: 

(XfXf ) = 5^ XrXfPMc(X„X,) . (49) 

Additionally, we define the variance {{AXiY) and the covariance (AXiAXj) respectively as: 

{{AXif) ^ (Xf) - (X,)2 , and (50) 
(AX, AX,) = (X,X,) - (X,)(X,) . (51) 

In the following, we use the scaled variance ooi and the correlation coefficient pij defined as: 

00, = ^^^^y^ , and (52) 
. , . (53) 



Let us consider a static and neutral system with four-velocity -u^ = (1,0,0,0), chemical 
potentials /i, = (0,0,0), local temperature T = = 0.160Gel^, and volume Vi = 2000/m'^. 
This is a system large enough^ for using the large volume approximation worked out in Sec- 
tion HH 

In Figs.([T]) and ([2]) we show the results of Monte Carlo runs of 2.5-10^ events each. Each value 
of A has been sampled 20 times to allow for calculation of a statistical uncertainty estimate. 
19 different values of A have been studied. In this case study, the extensive quantities baryon 
number B, strangeness S, electric charge Q, energy E, and longitudinal momentum Pz are 
considered for re-weighting. Conservation of transverse momenta Px and Py can be shown not 
to affect the ApT.i and Ai/i dependence of multiplicity fluctuations and correlations studied in 
the following sections. Their Atji dependence is, however, rather sensitive to P^ conservation. 



Angular correlations (not studiec 



joint Px and P„ conservation 



21 



in this article), on the other hand, are strongly sensitive to 



221. 



^ Generally it is not easy to say when a system is 'large enough' for the large volume approximation to be valid. 
Here we find good agreement with asymptotic analytic solutions. Charged systems, or Bose-Einstein/Fermi- 
Dirac systems, usually converge more slowly to their asymptotic solution. 
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FIG. 1: Mean values (left) and variances (right) of various extensive quantities, as listed in the legends, 
as a function of A. Each marker and its error bar represents the result of 20 Monte Carlo runs of 
2.5 • 10^ events each. 19 different equally spaced values of A have been investigated. Solid lines indicate 
GCE values (left), or linear extrapolations from the GCE value to the MCE limit (right). 
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FIG. 2: Covariances (left) and correlation coefficients (right) between various extensive quantities, as 
listed in the legends, as a function of A. Solid lines indicate linear extrapolations from the GCE value 
to the MCE limit (left), or GCE values (right). The rest as in Fig.©. 

In Fig.([T]) (left) we show the results for mean values of baryon number (B), strangeness (S), 
electric charge (Q), energy (E), and the momenta (P^) and (Pz). The solid lines represent 
GCE values. Only the expectation value of energy is not equal to 0, as the system sampled 
is assumed to be static and neutral with T ^ 0. The evolution of the respective variances 
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is shown in Fig.([T]) (right). Variances of extensive quantities subject to re- weighting converge 
hnearly to as A goes to 1. One notes that ((AP^.)^) remains constant (within error bars), as 
this quantity is not re-weighted in this case study. Please note that on many data points the 
error bars are smaller than the symbol used. 

In Fig.® (left) we show the evolution of covariances (ABAS), (ABAQ), (ASAQ), and 
{AEAQ) with the 'size of the bath'. As seen, the covariances between quantities considered 
for re- weighting also converge linearly to 0. In a neutral system, covariances between energy 
and charge are equal to 0. As an example, we show (AEAQ). In a static system, also the 
covariances between momenta and any other extensive quantity are equal to 0. As an example, 
we show (AEAPz). The correlation coefficients, Eq. fl53|) . on the other hand, remain constant 
as a function of A, as shown in Fig. ([2]) (right). The values of fully phase space integrated 
correlation coefficients pbs, Pbq, and psq can be compared to the GCE results denoted by the 
solid lines shown in Figs.(l5]- [7j) in Section IVll 

The variances and covariances converge linearly from their GCE values to their respective 
MCE limits in the large volume limit. The dependence of ((AXj)^), Eq. flSOl) . and (AXiAXj) , 
Eq. (l5Tl) . on the size of the bath A is given by: 

((AX,)2), = (1-A) ((AX,)2),,, + A ((AX,)2)_, (54) 
(AX.AX,), = (1 - A) {AX,AX,)g,, + A {AX,AX,)mce • (55) 

Mean values {Xi)x remain constant. This implies that the scaled variance u; of multiplicity 
fluctuations, Eq. fl52|) . also converges linearly: 

_ {(AN,)^), 

= Tl^M = (1 - ^) ^gce + A Umce , (56) 



from its GCE value cugce to the MCE limit Umce- Please note t 



lat Eqs. fl54|55|56l) are equivalent 



to the 'acceptance scaling' approximation^ used in 32|, |33|, |3J| . For the correlation coefficient, 

v/((AX,)2),((AX,)^)a ' 



4 



For the situation discussed here one could equivalently say that particles are randomly drawn from coordinate 



space of the total volume Vg. For the derivation of the acceptance scaling formula Sj] it was, however, assumed 
that particles are randomly drawn from a sample in momentum space. 
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the story is more complicated. In case both Xi and Xj are re-weighted and measured in full 
phase space, we find: 

{iAX,f)mce = {{AX^f)mce = {AX.AX,)^,^ = , (58) 

and the correlation coefficient px, Eg. (1571) . is independent of the value of A, see Fig.(l2j). In all 
other cases, one needs to extrapolate Eqs .( 1541551) separately, and then calculate the correlation 
coefficient. 

We have therefore successively transformed our Monte Carlo sample. As A ^ 1, we give 
larger and larger weight to events in the immediate vicinity of the equilibrium expectation value, 
and smaller and smaller weight to events away from it. The distribution of extensive quantities 
considered for re-weighting (a multivariate normal distribution in the GCE in the large volume 
limit) hence gets contracted to a 5- function with vanishing variances and covariances. I.e., 
we successively highlight the properties of events which have very similar values of extensive 
quantities. This will have a bearing on charge correlations and, in particular, multiplicity 
fluctuations and correlations discussed in the following sections. 
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FIG. 3: First and second moment of the weight factor Eq. pSI) as a function of A. The rest as in 
Fig.©. 

The price we pay is that, as A grows, so too does the statistical uncertainty. In the limit 
A ^ 1, we approach a sample-reject type of formalism. We cannot, therefore, directly obtain 
the microcanonical limit for the large system size studied here, as this is prohibited by available 
computing power. On the bright side, however, we can extrapolate to this limit. In Fig. ([3]) 
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we show the second moment of the weight factor, Eq. (l28l) . as a function of A. A large second 
moment (W^) imphes a large statistical uncertainty and, hence, usually requires a larger sample. 
We mention in this context that the intermediate ensembles, between the limits of GCE and 
MCE, may also be of phenomeno logical interest. 



V. MOMENTUM SPECTRA 



We next consider momentum spectra. In Fig.(jl]) we show transverse momentum and rapidity 
spectra of positively charged hadrons, both primordial and final state, for a static thermal 
system. 



Basec 
23| and 



on these momentum spectra we construct acceptance bins Ap^.i and A?/j, as in 21 



22 



35 



36l |. Momentum bins are constructed such that each of the five bins constructed 
contains on average one fifth of the total yield of positively charged particles. The values 
defining the bounds of the momentum space bins Ap^ j and A?/j are summarized in Table [H 




1.5 2 , , 2.5 

Pt [Gel/ 



FIG. 4: (Left:) Transverse momentum spectrum of positively charged hadrons, both primordial and 
final state. (Right:) Rapidity spectrum of positively charged hadrons, both primordial and final state. 
2 • 10^ events have been sampled. 

Resonance decay shifts the transverse momentum distribution to lower average transverse 
momentum (pj-) and widens the rapidity distribution of thermal 'fireballs' 



Final state 



transverse momentum bins are, hence, slightly 'contracted', while final state rapidity bins get 
slightly 'wider', when compared to their respective primordial counterparts. 
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PT,i [GeV] 


Pt,2 [GeV] 


PT,3 [GeV] 


PTA [GeV] 


PT,5 [GeV] 


PT,Q [GeV] 


primordial 


0.0 


0.22795 


0.36475 


0.51825 


0.73995 


5.0 


final state 


0.0 


0.17105 


0.27215 


0.38785 


0.56245 


5.0 




yi 


y2 




y4 


ys 


ye 


primordial 


-5.0 


-0.4275 


-0.1241 


0.1241 


0.4273 


5.0 


final state 


-5.0 


-0.5289 


-0.1553 


0.1551 


0.5289 


5.0 



TABLE I: Transverse momentum and rapidity bins Apx.i = [pT,i,PT,i+i] and Ayi = [yj,yj_|_i], both 
primordial and final state, for a static neutral Boltzmann system with temperature T = O.lGOGeV . 



Resonance decay combined with transverse as well as longitudinal flow is believed to provide 



a rather good description of experimental 
ion collisions at SPS and RHIC energies 



ly observed momentum spectra in relativistic heavy 



29 



30 



3l|. Our spectra, on the other hand, contain 



no flow and our results thus cannot be directly compared to experimental data or transport 
simulations. However, qualitatively one might observe effects of the kind discussed in the 
following. 



VI. THE MOMENTUM SPACE DEPENDENCE OF CORRELATIONS BETWEEN 
CONSERVED CHARGES 

An interesting example of quantities for which the measured value depends on the observed 
part of the momentum spectrum are the correlation coefficients between the charges baryon 
number B, strangeness S and electric charge Q. Please note that also the variances and 
covariances of the baryon number, strangeness, and electric charge distribution are sensitive 
to the acceptance cuts applied. Their values are additionally rather sensitive to the effects of 
globally enforced conservation laws. If the size of the 'bath' is reduced, a change in one interval 
of phase space will have to be balanced (preferably) by a change in another interval, and not 
by the 'bath'. 
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A. Grand Canonical Ensemble 

We will now consider the correlation coefficients pssi Pbqi and psq in limited acceptance 
bins Ap^ j and A|/j, as defined in Tabled in the grand canonical ensemble. Particles in one mo- 
mentum bin are then essentially sampled independently from particles in any other momentum 
space segment, due to the 'infinite bath' assumption. Nevertheless, the way in which quantum 
numbers are correlated is different in different momentum bins, as different particle species 
have, due to their different masses, different momentum spectra. 

Let us ffist make some basic observations about the hadron resonance gas and the way 
in which quantum numbers are correlated in a GCE. Charge fluctuations directly probe the 
degrees of freedom of a system, i.e. they are sensitive to its particle mass spectrum (and its 
quantum number configurations). We ffist consider the contribution of different particle species 
to the covariance {AXiAXj) , Eq.(l5Tl). and hence to the correlation coefficient pij, Eq.(l53l). 

All baryons have baryon number b = +1. Baryons can only carry strange quarks, i.e. their 
strangeness is always s < 0. Anti-baryons have b = —1, and s > 0. Hence, both groups 
contribute negatively to the baryon-strangeness covariance, and so {ABAS) < 0, and therefore 
P_B5 < 0, as indicated by the solid lines in Fig. ([5]). 

Positively charged baryons and their anti-particles contribute positively to the baryon- 
electric charge covariance (ABAQ) , while negatively charged baryons (and their anti-particles) 
contribute negatively. Two observations can be made on the hadron resonance gas mass spec- 
trum: there are more positively charged baryons than negatively charged ones, and their average 
mass is lower. I.e., in a neutral gas {pb = f^Q = Ps = ^) the contribution of positively charged 
baryons dominates and therefore (ABAQ) > and p^Q > 0, as indicated by the solid lines in 

Fig.m- 

Mesons and their anti-particles always contribute positively to the strangeness-electric charge 
correlation coefficient psq- Electrically charged strange mesons are either composed of a u- 
quark and an s-quark, or of an u-quark and a s-quark (and superpositions thereof). Their 
contribution to (ASAQ) is in either case positive. On the baryonic side, only the S'^ (as 
well as its degenerate states and their respective anti-particles) has a negative contribution to 
{ASAQ), while all other strangeness carrying baryons have either electric charge q = —1, or 
q = 0. Therefore, we find psq > 0, as indicated by the solid lines in Fig. ([7]). 

In Figs.([5llZl) ^6 show the correlation coefficients pss (baryon number - strangeness), Pbq 
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(baryon number - electric charge), and psq (strangeness - electric charge) as measured in the 
acceptance bins Ap^.i and A?/j defined in Table [U both primordial and final state. The average 
baryon number, strangeness, and electric charge in each bin is equal to zero, as the system 
is assumed to be neutral. The analytical primordial values (15 bins) shown in Figs.([5llZD 
calculated using analytical spectra. Please note that, again, on many data points the error bars 
are smaller than the symbol used. 
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FIG. 5: Baryon-strangeness correlation coefficient pBS in the GCE in limited acceptance windows, 
both primordial and final state. {Left:) transverse momentum bins Apj-.i- (Right:) rapidity bins Ayj. 
Horizontal error bars indicate the width and position of the momentum bins (And not an uncertainty!). 
Vertical error bars indicate the statistical uncertainty of 20 Monte Carlo runs of 10^ events each. The 
marker indicates the center of gravity of the corresponding bin. The solid lines show the fully phase 
space integrated GCE result. 

In Tables |TT] to [IV] we summarize the transverse momentum and rapidity dependence of 
the correlation coefficients pbs, Pbq, and psQ- The statistical error quoted corresponds to 20 
Monte Carlo runs of 10^ events each. The analytical values (5 bins) listed in the tables are 
calculated using the momentum bins defined in Table [H Mild differences between Monte Carlo 
and analytical results are unavoidable. The analytical values are also not exactly symmetric 
in Ayi, as the exact size of the acceptance bins constructed is sensitive to the number of bins 
used for the calculation of the momentum spectra. The values of the correlation coefficient p 
are also rather sensitive to exact bin size, and the fourth digit becomes somewhat unreliable. 

We next attempt to explain, in turn, the rapidity dependence of pss, Pbq, and psg- Strange 
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FIG. 6: Baryon-electric charge correlation coefficient pBQ in the GCE in hmited acceptance windows, 
both primordial and final state. (Left:) transverse momentum bins Apx.i- (Right:) rapidity bins Ayi. 
The rest as in Fig. ([5]). 
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FIG. 7: Strangeness-electric charge correlation coefficient psQ in the GCE in limited acceptance 
windows, both primordial and final state. (Left:) transverse momentum bins ApT,i- (Right:) rapidity 
bins Ayi. The rest as in Fig. ([5]). 

baryons are, on average, heavier than non-strange baryons, so their rapidity distributions 
are narrower. The kaon rapidity distribution is then, compared to baryons, again wider. A 
change in baryon number (strangeness) at high |y| is less hkely to be accompanied by a change 
in strangeness (baryon number) than at low \y\. The value of Pbs, therefore, drops toward 
higher rapidity, as shown in Fig. ([5]), (right). By the same argument, we find a weakening 
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PBS 




ApT,2 


ApT,3 


ApT,4 


ApT,5 


calc 

rpVlTfl 


-0.2479 


-0.2641 


-0.2864 


-0.3188 


-0.3839 


Pprim 


-0.248 ± 0.003 


-0.264 ± 0.003 


-0.286 ± 0.003 


-0.319 ±0.002 


-0.385 ± 0.002 


P final 


-0.216 zb 0.002 


-0.220 ±0.003 


-0.241 ± 0.004 


-0.269 ±0.003 


-0.335 ± 0.003 


PBS 




Ay2 


Ays 


Ay4 


Ays 


„calc 
Pprim 


-0.2407 


-0.3345 


-0.3536 


-0.3345 


-0.2408 


Pprim 


-0.241 ± 0.003 


-0.334 ± 0.003 


-0.353 ±0.003 


-0.335 ± 0.003 


-0.240 ± 0.003 


P final 


-0.191 ±0.002 


-0.300 ± 0.002 


-0.328 ± 0.002 


-0.299 ±0.002 


-0.190 ±0.002 



TABLE II: Baryon-strangeness correlation coefficient pss in the GCE in transverse momentum bins 
Apr.i and rapidity bins Ay^, both primordial and final state. For comparison, analytical values /O^"'^ 
for primordial correlations are included. The statistical uncertainty corresponds to 20 Monte Carlo 
runs of 10^ events each. 



PBQ 


ApT,l 


ApT,2 


ApT,3 


ApT,4 


ApT,5 


^calc 
Pprim 


0.1120 


0.1271 


0.1420 


0.1579 


0.1781 


Pprim 


0.113 ±0.002 


0.126 ± 0.002 


0.143 ±0.003 


0.158 ±0.002 


0.178 ± 0.003 


P final 


0.112 ±0.003 


0.120 ±0.003 


0.138 ±0.003 


0.164 ±0.003 


0.221 ±0.003 


PBQ 


Ayi 


Ay2 


Ays 


Ay4 


Ays 


^calc 
Pprim 


0.1160 


0.1601 


0.1658 


0.1601 


0.1160 


Pprim 


0.116 ±0.002 


0.160 ±0.003 


0.166 ±0.003 


0.159 ±0.003 


0.117 ±0.002 


P final 


0.118 ±0.003 


0.192 ± 0.003 


0.202 ± 0.003 


0.192 ±0.003 


0.119 ±0.003 



TABLE III: Baryon-electric charge correlation coefficient pbq in the GCE in transverse momentum 
bins ApT^i and rapidity bins Ayj, both primordial and final state. 



of the baryon-electric charge correlation pbq at higher rapidity (Fig. ([6]), [right]) as the 
rapidity distribution of electrically charged particles is wider than that of baryons. For the 
strangeness-electric charge correlation coefficient we find first a mild rise, and then a somewhat 
stronger drop of psq towards higher rapidity. As one shifts ones acceptance window towards 
higher values of \y\, first the contribution of baryons (in particular E"*") decreases and, as the 
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PSQ 




ApT,2 


ApT,3 


ApT,4 


ApT,5 


calc 
rprim 


0.2831 


0.3033 


0.3150 


0.3185 


0.3055 


Pprim 


0.284 ± 0.003 


0.304 ± 0.003 


0.314 ±0.003 


0.319 ±0.002 


0.305 ± 0.002 


P final 


0.243 ± 0.003 


0.254 ± 0.003 


0.276 ± 0.003 


0.292 ±0.003 


0.303 ± 0.002 


PSQ 




Ay2 


Ays 


Ay4 


Ay5 


^calc 
Pprim 


0.2934 


0.3137 


0.3104 


0.3137 


0.2934 


Pprim 


0.294 ± 0.003 


0.314 ± 0.003 


0.310 ±0.002 


0.312 ±0.003 


0.292 ± 0.002 


P final 


0.255 ±0.002 


0.299 ± 0.003 


0.297 ±0.003 


0.298 ± 0.003 


0.255 ± 0.003 



TABLE IV: Strangeness-electric charge correlation coefficient psq in the GCE in transverse momen- 
tum bins ApT^i and rapidity bins Ayi, both primordial and final state. 

meson contribution grows, psQ rises slightly. Towards the highest \y\, pions again dominate 
and de-correlate the quantum numbers. 

The transverse momentum dependence can be understood as follows: heavier particles 
have higher average transverse momentum (px) and, hence, their influence increases towards 
higher pt- Heavy particles have a tendency to carry several charges, causing the correlation 
coefficients to grow. 

The contribution of strange baryons compared to non-strange baryons grows towards higher 
transverse momentum, as strange baryons have on average larger mass than non-strange 
baryons. The correlation coefficient pBs thus becomes strongly negative at high pt- As the 
contribution of baryons compared to mesons grows stronger towards larger pt, a change in 
baryon number (electric charge) is now more likely to be accompanied by a change in electric 
charge (baryon number) than at low Pt, and Pbq increases with px (The A resonances^ ensure 
it keeps rising). For the Ap^j dependence of psq we finally note that one of the strongest 
contributors at higher pj- is the Q~ , with a relatively low mass of rriQ- = 1.672GeV . So after 
a rise, psq drops again towards highest px, due to an increasing S"*" contribution^. 



^ Included in the THERMITS particle table up to the A (2420) . 
^ Included in the THERMUS particle table up to the E(2030). 
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Since resonance decay has the habit of dropping the hghter particles (mesons) at low px 
and higher \y\, while keeping heavier particles (baryons) at higher pt and at mid-rapidity, 
none of the above arguments about the transverse momentum and rapidity dependence are 
essentially changed by resonance decay. The correlation coefficient pbs becomes more negative 
towards higher px-, while becoming weaker towards higher \y\. Similarly, pbq grows larger at 
high Pt and drops towards higher y. The larger contributions of baryons to the high pt tail 
of the transverse momentum spectrum, and their decreased contribution to the tails of the 
rapidity distribution, compared to mesons, are to blame. The bump in the px dependence of 
PsQi presumably caused by the S"*", has vanished, as the S+ is only considered as stable in its 
lightest version with mass ms+ = l.lSOGeV. The small bump in the y dependence of psq-, 
however, stays. The correlation is presumably first increased by a growing kaon contribution 
and then again decreased by a growing pion contribution at larger rapidities. 

The values of p after resonance decay are directly sensitive to how the data is analyzed. 
In the above study we analyzed final state particles (stable against strong decays) only. One 
could, however, also reconstruct decay positions and momenta of parent resonances and could 
then count them as belonging to the acceptance bin the parent momentum would fall into. In 
the situation above, however, this would again yield the primordial scenario. If reconstruction 
of resonances is not done, one is sensitive to charge correlations carried by final state particles. 
As in the primordial case, a larger acceptance bin effectively averages over smaller bins. How- 
ever, the smaller the acceptance bin, the more information is lost due to resonance decay. In 
full acceptance, final state and primordial correlation coefficients ought to be the same, since 
quantum numbers (and energy-momentum) are conserved in the decays of resonances. 

B. Extrapolating to the MCE 

We next consider the extrapolation to the MCE limit of variances and covariances and, hence, 
correlation coefficients, of joint distributions of charges in limited acceptance. The primordial 
joint baryon number - strangeness distributions in different transverse momentum bins will 
serve as examples. In this subsection, we use an extended data set of 20 ■ 8 ■ 10^ events. 

In Fig. ([8]) we show the evolution of the variances of the marginal primordial baryon number 
distribution ((Ai?)^) {left) and of the marginal primordial strangeness distribution ((AS*)^) 
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FIG. 8: Evolution of the variance of the marginal baryon number distribution {{AB)'^) (left) and 
the variance of the marginal strangeness distribution ((AS")^) (right) with A for a primordial hadron 
resonance gas in different ApT.i bins. Each marker and its error bar except the last represents the result 
of 20 Monte Carlo runs of 10^ events each. 8 different equally spaced values of A have been investigated. 
The last marker denotes the result of the extrapolation. Solid lines indicate extrapolations from the 
GCE value to the MCE hmit. 
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FIG. 9: Evolution of the covariance (ABAS) (left) and the correlation coefficient pss {fight) of the 
baryon number - strangeness distribution with A for a primordial hadron resonance gas in different 
ApT,i bins. The rest as in Fig. ([8]). 



[right) in the transverse momentum bins Ap^- j, defined in Table [H as a function of the size 
of the bath X = Vi/Vg. 8 equally spaced values of A have been investigated. The last marker 
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denotes the result of the extrapolation. In Fig. ([9]) we show the dependence of the primordial 
covariance (ABAS) (left) and the primordial correlation coefficient pbs (right) of the joint 
baryon number - strangeness distribution on the size of the bath A. 

Let us first comment on the GCE values of variances (the left most markers in Fig.®). 
As each of the 5 momentum bins holds one fifth of the charged particle yield and, hence, less 
than one fifth of the baryonic contribution in the lowest bin Apx.i, and more than one fifth 
in the highest bin ApT,5, we find the baryon number variance ((Ai?)^) largest in ApT,5, and 
smallest in Ap-r^i- If binned in rapidity: Ay^ has the strongest baryon contribution, and, hence, 
{{AB)"^) is largest there. The same goes for the variance {{AS)^) of the marginal strangeness 
distribution. Strangeness carrying particles are on average heavier than electrically charged 
particles and, hence, the strangeness contribution is strongest around mid-rapidity and towards 
larger transverse momentum (i.e. ((AS*)^) is largest in Ay^ and Ap^.s, while being smallest in 
Ay I, Ay 5, and ApT,i)- 
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FIG. 10: MCE baryon number - strangeness correlation coefficient pBS in limited acceptance windows, 
both primordial and final state. (Left:) transverse momentum bins Apx.i- (Right:) rapidity bins Ay^. 
Horizontal error bars indicate the width and position of the momentum bins (And not an uncertainty!). 
Vertical error bars indicate the statistical uncertainty of the extrapolation of 8 • 20 Monte Carlo runs 
of 10^ events each. The marker indicates the center of gravity of the corresponding bin. The solid 
lines show the fully phase space integrated GCE result. 



The ApT j dependence of the GCE covariance (Ai^AS*) and the GCE correlation coefficient 
Pbs in Fig. ([9]) is explained by the arguments of the previous subsection. Varying contributions 
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of hadrons of different mass (and cliarge contents) to different parts of momentum space are 
responsible. 

We now turn our attention to the extrapolation. MCE effects on the baryonic sector are felt 
most strongly in momentum space segments in which the baryonic contribution is strong (e.g. 
see the evolution of the last bin Apj-^^ with A in Figs. flHHQ]) ). The correlation coefficient is not 
as strongly affected, in general, by MCE effects. 

In Fig. ( fTOl) we show the results of the extrapolation to the MCE limit of the baryon number- 
strangeness correlation coefficient pss in acceptance bins Apx^i and Ayi, both primordial and 
final state. MCE values are closer to each other than corresponding GCE values. Fig. ([5]). The 
influence of globally applied conservation laws on charge correlations is less strong than for the 
multiplicity fluctuations and correlations discussed in the next section. 

VII. MOMENTUM SPACE DEPENDENCE OF MULTIPLICITY FLUCTUATIONS 
AND CORRELATIONS 

Multiplicity fluctuations and correlations are qualitatively affected by the choice of ensemble 
and are directly sensitive to the fraction of the system observed. For vanishing size of ones 
acceptance window, one would lose all information on how the multiplicities of any two distinct 
groups Ni and Nj of particles are correlated, and measure = 0. This information, on the 
other hand, is to some extent preserved in pss, Pbq, and psq, i.e. the way in which quantum 
numbers are correlated, if at least occasionally a particle is detected during an experiment. 

We first sample the same GCE system, which we have discussed in the previous sections, 
and consider the effects of resonance decay. Next the joint distributions of positively and 
negatively charged particles in momentum bins Apj-^i and Arji are constructed. Then we, in 
turn, extrapolate the GCE primordial and final state results on the scaled variance u, Eq. (!52|) . 
and the correlation coefficient p, Eq. (!53l) . to the MCE limit. 

A. Grand Canonical Ensemble 

In Fig.f fTTl) we show the ApT,i {left) and Ayi (right) dependence of the GCE scaled variance 
UJ+ of positively charged hadrons, both primordial and final state. In the primordial Boltz- 
mann case one finds no dependence of multiplicity fluctuations on the position and size of the 
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acceptance window. The observed multiplicity distribution is, within error bars, a Poissonian 
with scaled variance a;+ = 1. In fact, in the primordial GCE Boltzmann case any selection of 
particles has u = 1. 
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FIG. 11: GCE scaled variance 0;+ of multiplicity fluctuations of positively charged hadrons, both 
primordial and final state, in transverse momentum bins Apx^i {left) and rapidity bins Ayi (right). 
Horizontal error bars indicate the width and position of the momentum bins (And not an uncertainty!). 
Vertical error bars indicate the statistical uncertainty of 20 Monte Carlo runs of 2 • 10^ events each. 
The markers indicate the center of gravity of the corresponding bin. The solid line indicates the final 
state acceptance scaling estimate. 

In Fig. (fT2l) we show the Ap^^i (left) and Ayi (right) dependence of the GCE correlation 
coefficient p_|__ between positively and negatively charged hadrons, both primordial and final 
state. In the primordial Boltzmann case one finds also no dependence of multiplicity correlations 
on the position and size of the acceptance window. The observed joint multiplicity distribution 
is a product of two Poissonians with correlation coefficient p_) = 0. 

Resonance decay is the only source of correlation in an ideal GCE Boltzmann gas. Neutral 
hadrons decaying into two hadrons of opposite electric charge are the strongest contributors to 
the correlation coefficient p+-. The chance that both (oppositely charged) decay products are 
dropped into the same momentum space bin is obviously highest at low transverse momentum 
(i.e. the correlation coefficient is strongest in ApT,i)- The rapidity dependence is somewhat 
milder again, because heavier particles (parents) are dominantly produced at mid-rapidity and 
spread their daughter particles over a range in rapidity. One notes that the scaled variances 
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FIG. 12: GCE multiplicity correlations p-\ between positively and negatively charged hadrons, both 

primordial and final state, in transverse momentum bins Apx^i {left) and rapidity bins Ayi (right). 
The rest as in Fig. (fTT]) . 

and correlation coefficients in the respective acceptance bins in Figs.( fTT][T2l) are generally larger 
than the acceptance scaling procedure^ suggests, with the notable exception of p+^{ApT,5). 

If one would construct now a larger and larger number of momentum space bins of equal 
average particle multiplicities, one would successively lose more and more information about 
how multiplicities of distinct groups of particles are correlated. 

There is a simple relation connecting the scaled variance of the fluctuations of all charged 
hadrons uj± to the fluctuations of only positively charged particles via the correlation coef- 
ficient p_|__ between positively and negatively charged hadrons in a neutral system: 

= ^+ (1 + p+_). (59) 

We, therefore, find the effect of resonance decay on the Apx^i dependence of uj±. to be consid- 
erably stronger than on that of u;+, and generally uj± > w+, as the correlation coefficient p+_ 
remains positive in the final state GCE. Compared to this, the final state values of u!±, and 
remain rather fiat with Ai/i in the GCE. 

^ For the acceptance scaling approximation it is assumed that particles are randomly detected with a certain 
probability q = 0.2, independent of their momentum. 
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B. Extrapolating to the MCE 



In the very same way that we extrapolated fully phase space integrated extensive quantities 
to the MCE limit in Section [TVl we now extrapolate multiplicity fluctuations uoj^ and correlations 

in transverse momentum bins Ap-p j and rapidity bins Ai/j for a hadron resonance gas from 

the GCE (A = 0) to the MCE fA — > 1). Analytical primordial MCE results are done in the 



infinite volume approximation [2l|, |22]. We, hence, have some guidance as to further asses the 
accuracy of the extrapolation scheme. For final state fluctuations and correlations in limited 
acceptance, on the other hand, no analytical results are available. 

Mean values of particle numbers of positively charged hadrons (A^+) and negatively charged 
hadrons (iV„) in the respective acceptance bins, defined in Tabled remain constant as A goes 
from to 1, while the variances ((AA^+)^) and ((AA^„)^), and covariance (AA^+AA^-) converge 
linearly to their respective MCE limits. The correlation coefficient p+_ between positively and 
negatively charged hadrons, on the other hand, will not approach its MCE value linearly, as 
discussed in Section HVl 

1. Primordial 

In Fig. (fT3|) we show the primordial scaled variance u;+ of positively charged hadrons in 
transverse momentum bins Ap^ j {left) and rapidity bins A?/j [right) as a function of the size 
of the bath A = Vi/Vg, while in Fig. (fT4|) we show the dependence of the primordial correlation 
coefficient p+_ between positively and negatively charged hadrons in transverse momentum 
bins Apr.j [left) and rapidity bins Ay^ [right] on A. 

The results of 8 ■ 20 Monte Carlo runs of 2 ■ 10^ events each are summarized in Table IVl 
The system sampled was assumed to be neutral pj = (0,0,0) and static = (1,0,0,0) with 
local temperature (3~^ = O.lGOGeV^ and a system volume of Vi = 2000 fm^. 8 different values 
of A have been studied. The last marker (A = 1) denotes the result of the extrapolation. Only 
primordial hadrons are analyzed. Values for both Apx.i and Avi bins are listed. Analytical 



numbers are calculated according to the method developed in [2l|, l22|, using the acceptance 
bins defined in Table HI and are shown for comparison. 

The effects of energy-momentum and charge conservation on primordial multiplicity fluctua- 



tions and correlations in finite acceptance have been discussed in 
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22|. A few words attempt 
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to summarize. 
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FIG. 13: Evolution of the primordial scaled variance ujj^ of positively charged hadrons with the Monte 
Carlo parameter A = Vi/Vg for transverse momentum bins ApT,i (left) and for rapidity bins Am 
(right). The solid lines show an analytic extrapolation from GCE results (A = 0) to the MCE limit 
(A ^ 1). Each marker and its error bar except the last represents the result of 20 Monte Carlo runs of 
2 • 10^ events. 8 different equally spaced values of A have been investigated. The last marker denotes 
the result of the extrapolation. 
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FIG. 14: Evolution of the primordial correlation coefficient between positively and negatively 

charged hadrons with the Monte Carlo parameter A = Vi/Vg for transverse momentum bins Apx^i 
(left) and for rapidity bins Ayi (right). The rest as in Fig. (fT3]) . 
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Let us first attend to fully phase space integrated results. The scaled variance of multi- 
plicity fluctuations is lowest in the MCE due to the requirement of exact energy and charge 
conservation, somewhat larger in the CE, and largest in the GCE, as now all constraints on 



the microstates of the system have been dropped [27|, l33|, l3J] . The fully phase space integrated 
MCE and CE correlation coefficients between oppositely charged particles are rather close to 
1. Doubly charged particles allow for mild deviation, as also the A^"*" resonance is counted as 
only one particle. 

The transverse momentum dependence can be understood as follows: a change in particle 
number at high transverse momentum involves a large amount of energy. I.e., in order to balance 
the energy record, one needs to create (or annihilate) either a lighter particle with more kinetic 
energy, or two particles at lower pt- This leads to suppressed multiplicity fluctuations in high 
ApT,i bins compared to low ApT,i bins. By the same argument, it seems favorable, due to 
the constraint of energy and charge conservation, to balance electric charge, by creating (or 
annihilating) pairs of oppositely charged particles, predominantly in lower Ap^ j bins, while 
allowing for a more un-correlated multiplicity distribution, i.e. also larger net-charge {6Q = 
— N^) fluctuations, in higher Ap^^i bins. 

For the rapidity dependence similar arguments hold. Here, however, the strongest role 
is played by longitudinal momentum conservation. A change in particle number at high y 
involves now, in addition to a large amount of energy, a large momentum pz to be balanced. 
The constraints of global conservation are, hence, felt least severely around \y\ ~ 0, and it 
becomes favorable to balance charge predominantly at mid-rapidity (p+_ larger) and allow for 
stronger multiplicity fluctuations (u;_|_ larger) compared to forward and backward rapidity bins. 

In a somewhat casual way one could say: events of a neutral hadron resonance gas with 
values of extensive quantities B, S, Q, E and Pz in the vicinity of {Q{) have a tendency to have 
similar numbers of positively and negatively charged particles at low transverse momentum pt 
and rapidity y and less strongly so at high px and \y\. 

The statistical error on the 'data' points grows as A ^ 1, as can be seen from Figs. fll3lll4p . 
The extrapolation helps greatly to keep the statistical uncertainty on the MCE limit low, as 
summarized in Table |Vl and can be seen from a comparison of the last two data points in 
Figs.l ITSfT^ . The last point and its error bar denote the result of a linear extrapolation of 
variances and covariances, while the second to last data point and its error bar are the result 
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primordial 


A??T 1 


Ax>T 9 




Al*T^ A 


— 'x-^-t ,0 




1.000 ± 0.002 


1.000 ± 0.002 


1.000 ± 0.002 


1.000 ± 0.002 


1.000 ± 0.002 


, ,mce 


0.889 ± 0.007 


0.880 ± 0.007 


0.869 ± 0.007 


0.850 ± 0.006 


0.798 ± 0.007 


mce,c 


0.8886 


0.8802 


0.8682 


0.8489 


0.7980 




0.000 ± 0.002 


-0.000 ± 0.002 


-0.000 ± 0.002 


0.000 ± 0.002 


0.000 ± 0.001 




0.094 ± 0.005 


0.085 ± 0.006 


0.072 ± 0.006 


0.056 ± 0.006 


0.003 ± 0.005 


mce.c 
1 -] — 


0.0935 


0.0844 


0.0730 


0.0554 


0.0040 


primordial 






ATy^ 




Ay^ 


^+ 


1 000 ± 002 


1 000 ± 002 


1 000 ± 003 


1 000 ± 002 


1 000 ± 002 


nice 


0.795 ± 0.006 


0.835 ± 0.007 


0.853 ± 0.008 


0.834 it 0.006 


0.794 ± 0.007 


mce,c 


0.7950 


0.8350 


0.8521 


0.8351 


0.7949 


nee 
P+- 


-0.000 ± 0.001 


0.000 ± 0.002 


0.001 ± 0.002 


0.000 ± 0.002 


-0.000 ± 0.002 


„mce 
P+- 


-0.013 ±0.005 


0.040 ± 0.006 


0.061 ± 0.006 


0.041 ± 0.006 


-0.012 ±0.006 


mee,c 
P+~ 


-0.0135 


0.0406 


0.0616 


0.0406 


-0.0135 



TABLE V: Summary of the primordial scaled variance a;+ of positively charged hadrons and the 

correlation coefficient between positively and negatively charged hadrons in transverse momentum 

bins i^PT,i and rapidity bins Ay,,. Both the GCE result (A = 0) and the extrapolation to MCE (A = 1) 
are shown. The uncertainty quoted corresponds to 20 Monte Carlo runs of 2 • 10^ events (GCE) or is 
the result of the extrapolation (MCE). Analytic MCE results uj^^^'^ and p^^"^'^ are listed too. 



of 20 Monte Carlo runs with A = 0.875. The analytical MCE values are well within error bars 
of extrapolated Monte Carlo results, and agree surprisingly well, given the large number of 
"conserved" quantities (5) and a relatively small sample size of 8 ■ 20 ■ 2 ■ 10^ = 3.2 ■ 10'' events. 
In a sample-reject type of approach this sample size would yield a substantially larger statistical 
error, as only events with exact values of extensive quantities are kept for the analysis. As the 
system size is increased, a sample-reject formalism, hence, becomes increasingly inefficient, while 
the extrapolation method still yields good results. For a further discussion see Appendix \^ 
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2. Final State 



We now attend to the extrapolation of final state multiplicity fluctuations and correlations 
to the MCE limit. An independent Monte Carlo run for the same physical system was done, 
but now with only stable final state particles 'detected'. 

In Fig.l fTSl) we show the final state scaled variance uj+ of positively charged hadrons in 
transverse momentum bins ^pT,i (left) and rapidity bins Ayi (right) as a function of A, while in 
Fig. (fT6l) we show the dependence of the final state correlation coefficient p+_ between positively 
and negatively charged hadrons in transverse momentum bins ApT^i (left) and rapidity bins Ai/i 
(right) on the size of the bath A = Vi/Vg. 
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A A 

FIG. 15: Evolution of the final state scaled variance uj^ of positively charged hadrons with the Monte 
Carlo parameter A = Vi/Vg for transverse momentum bins Apj^i (left) and for rapidity bins Ayi 
(right). The solid lines show an analytic extrapolation from GCE results (A = 0) to the MCE limit 
(A — > 1). Each marker except the last represents the result of 20 Monte Carlo runs of 2 • 10^ events. 
8 different equally spaced values of A have been investigated. The last marker denotes the result of 
the extrapolation. 

The ApT^i and Ayi dependence on A of the final state MCE scaled variance is qualitatively 
similar to that of the primordial versions, Fig. flT^ . and is essentially also explained by the 
arguments of the previous section. The effects of charge and energy-momentum conservation 
work in pretty much the same way as before, and it still seems favorable to have events with 
wider multiplicity distributions at low pj- and low y, and narrower distributions at larger p^ 
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FIG. 16: Evolution of the final state correlation coefficient pj^ between positively and negatively 

charged hadrons with the Monte Carlo parameter A = Vi/Vg for transverse momentum bins ApT^i 
(left) and for rapidity bins Ayi (right). The rest as in Fig. (|15p . 



and larger \y\. The dependence of the final state correlation coefficients p+_ on A, Fig. f|T6|) . is 
a bit different to the primordial case, Fig. (1141) . However, in the MCE limit, events still tend to 
have more similar numbers of oppositely charged particles at low pt and low y, than at large 
Pt and large \y\. 

The effects of resonance decay are qualitatively different in the MCE, CE, and GCE. Let us 



again first attend to fully phase space integrated multiplicity fiuctuations discussed in |33l. 134 1. 
The final state scaled variance increases in the GCE and CE compared to the primordial scaled 
variance. Multiplicity fiuctuations of neutral mesons remain unconstrained by conservation 
laws. However, they often decay into oppositely charged particles, which increases multiplicity 
fiuctuations of pions, for instance. In the MCE, due to the constraint of energy conservation, the 
event-by-event fiuctuations of primordial pions are correlated to the event-by-event fiuctuations 



of, in general, primordial parent particles, and uj^^'^"'^ < 



prim 



is possible in the MCE. 



In Fig. f[T7|) and Fig. flTS]) we compare the final state Ap^i (left) and Ayi (right) dependence 
of the MCE scaled variance and the MCE correlation coefficient p+_ respectively to their 
primordial counterparts. The results of 8 ■ 20 Monte Carlo runs of 2 ■ 10^ events each for a static 
and neutral hadron resonance gas with T = O.lGOGeV are summarized in Table fIVip . 

A few words to summarize Figs. (ll7fl8p : resonance decay and (energy) conservation laws 
work in the same direction, as far as the transverse momentum dependence of the scaled 
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FIG. 17: MCE scaled variance uj+ of multiplicity fluctuations of positively charged hadrons, both 
primordial and final state, in transverse momentum bins ApT^i (left) and rapidity bins Ayi (right). 
Horizontal error bars indicate the width and position of the momentum bins (And not an uncertainty!). 
Vertical error bars indicate the statistical uncertainty quoted in Table IVTl The markers indicate the 
center of gravity of the corresponding bin. The solid and the dashed lines show final state and 
primordial acceptance scaling estimates respectively. 
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FIG. 18: MCE multiplicity correlation coefficient between positively and negatively charged 

hadrons, both primordial and final state, in transverse momentum bins Apx.i (left) and rapidity bins 
Ayi (right). The rest as in Fig.(|17p. 



variance u;+ and the correlation coefficient p+_ is concerned. Both effects lead to increased 
multiplicity fluctuations and an increased correlation between the multiplicities of oppositely 
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final state 






APT,3 


ApT,4 


ApT,5 




1.031 ±0.002 


1.026 ±0.002 


1.020 ± 0.002 


1.015 ±0.002 


1.010 ±0.002 


, ,mce 


0.904 ± 0.007 


0.884 ± 0.007 


0.872 ± 0.007 


0.847 ±0.007 


0.778 ± 0.006 




0.163 ±0.001 


0.107 ±0.001 


0.109 ±0.001 


0.075 ± 0.002 


0.052 ± 0.002 


— 


0.143 ± 0.005 


0.088 ± 0.005 


0.090 ± 0.005 


0.049 ± 0.006 


-0.010 ±0.006 


final state 






Ays 


Ay4 


Ays 




1.017 ±0.002 


1.023 ±0.002 


1.024 ±0.002 


1.023 ±0.003 


1.017 ±0.002 


, mice 


0.771 ± 0.007 


0.840 ± 0.006 


0.859 ± 0.007 


0.839 ± 0.007 


0.770 ± 0.006 


pT- 


0.100 ±0.001 


0.116 ±0.001 


0.115 ±0.002 


0.115 ±0.002 


0.100 ±0.001 


P+- 


-0.027 ±0.005 


0.069 ± 0.005 


0.092 ± 0.006 


0.069 ± 0.006 


-0.027 ±0.005 



TABLE VI: Summary of the final state scaled variance w+ of positively charged hadrons and the 

correlation coefficient p_j between positively and negatively charged hadrons in transverse momentum 

bins ApT^i and rapidity bins Ay^. Both the GCE result (A = 0) and the extrapolation to the MCE 
(A = 1) are shown. The uncertainty quoted corresponds to 20 Monte Carlo runs of 2 • 10^ events 
(GCE) or is the result of the extrapolation (MCE). 



charged particles in the low px region, compared to the high pr domain. 

Compared to this, the MCE Ayi dependence of uj_^ and p+_ is mainly dominated by global 
conservation of P^. Resonance decay effects, see Figs. lfTTfT^ . are more equal across rapidity, 
than in transverse momentum. 

Again, we find the scaled variance of all charged particles larger than the scaled variance of 
only positively charged hadrons uj± > except for when < 0, i.e when the multiplicities 
of oppositely charged particles are anti-correlated, as for instance in ApT,5, Ai/i, and Ay 5. In 
contrast to that, we narrowly find cu-i- > 1 in the lowest transverse momentum bin Apx^i- 
The qualitative picture presented in Fig. (fT7|) could be compared to similar analysis of 



UrQMD transport simulation data [35 



tuations in limited momentum bins 



or recently published NA49 data on multiplicity fluc- 



3611. We, however, do not claim that the effects discussed 



above are the sole effects leading to the qualitative agreement with either of the two. 
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VIII. SUMMARY 



We have presented a recipe for a thermal model Monte Carlo event generator capable of 
extrapolating fluctuation and correlation observables for Boltzmann systems of large volume 
from their GCE values to the MCE limit. Our approach has a strong advantage compared to 
analytical approaches or standard microcanonical sample- and- reject Monte Carlo techniques, 
in that it can handle resonance decays as well as (very) large system sizes at the same time. 

To introduce our scheme, we have conceptually divided a microcanonical system into two 
subsystems. These subsystems are assumed to be in equilibrium with each other, and subject to 
the constraints of joint energy-momentum and charge conservation. Particles are only measured 
in one subsystem, while the second subsystem provides a thermodynamic bath. By keeping the 
size of the first subsystem fixed, while varying the size of the second, one can thus study the 
dependence of statistical properties of an ensemble on the fraction of the system observed (i.e. 
assess their sensitivity to globally applied conservation laws). The ensembles generated are 
thermodynamically equivalent in the sense that mean values in the observed subsystem remain 
unchanged when the size of the bath is varied, provided the combined system is sufficiently 
large. 

The Monte Carlo process can be divided into four steps. In the first two steps primordial 
particle multiplicities for each species, and momenta for each particle, are generated for each 
event by sampling the grand canonical partition function. In the third step resonance decay of 
unstable particles is performed. Lastly the values of extensive quantities are calculated for each 
event and a corresponding weight factor is assigned. All events with the same set of extensive 
quantities hence still have 'a priori equal probabilities'. In the limit of an infinite bath, all 
events have a weight equal to unity. In the opposite limit of a vanishing bath, only events 
with an exactly specified set of extensive quantities have non- vanishing weight. In between, 
we extrapolate in a controlled manner. The method is even rather efficient for large volume, 
inaccessible to sample-and-reject procedures, and agrees well, where available, with analytic 
asymptotic microcanonical solutions. 

Given the success of the hadron resonance gas model in describing experimentally measured 
average hadron yields, and its ability to reproduce low temperature lattice susceptibilities, the 
question arises as to whether fiuctuation and correlation observables also follow its main line. In 
particular, three effects are nicely discussable: Resonance decay, conservation laws, and limited 
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acceptance effects. Due to tlie Monte Carlo nature, data can be analyzed in close relation to 
experimental analysis techniques. The hadron resonance gas is an ideal testbed for this type of 
study, in that it is simple and intuitive. 

The statistical properties of a sample of hadron resonance gas events show a systematic 
dependence on what part of the momentum distribution and what fraction of the system is 
observed. Two examples served to illustrate: grand canonical charge-charge correlations, and 
microcanonical multiplicity fluctuations and correlations. In the case of charge-charge corre- 
lations, momentum space effects are caused by different masses of hadrons and, hence, their 
varying contribution to different parts of the momentum spectra. Although microcanonical 
effects on the (co)variances of the joint baryon number - strangeness - electric charge distribu- 
tion are considerable, they remain weak for the correlation coefficients between these quantum 
numbers. In contrast to this, momentum space effects on multiplicity fluctuations and correla- 
tions arise due to conservations laws. For an ideal primordial grand canonical ensemble in the 
Boltzmann approximation (our starting point), multiplicity distributions are just uncorrelated 
Poissonians, regardless of the acceptance cuts applied, as particles are assumed to be produced 
independently. The requirement of energy-momentum and charge conservation leads to sup- 
pressed fluctuations and enhanced correlations between the multiplicities of two distinct groups 
of particles at the 'high momentum' end of the momentum spectrum, provided some fraction of 
an isolated system is observed. Resonance decay does not change these trends. The arguments 
on which the explanation of this particular dependence are based seem general enough to hope 
that they might hold too in non-equilibrium systems, such as real heavy ion data or theoretical 
transport simulations. 

A direct comparison with experimental data seems problematic at the moment. The static 
global thermal and chemical equilibrium assumption made here is certainly insufficient. The 
model presented here is far from complete. Several interesting aspects deserve attention. They 
include the sampling of Fermi-Dirac or Bose-Einstein particles, for which low transverse mo- 
mentum is particularly sensitive; finite volume corrections could be done (possible if one has a 
good approximation to W); the convergence properties (at fixed A, and as a function of A) fall 
basically into the same direction; so far we also have not derived a thermodynamic potential for 
our ensembles; one could also consider more general forms of W; one could ask how to couple 
two systems of different densities, or altogether depart from the local equilibrium assumption. 



39 



There are also several interesting things that the model could do in its present form. Examples 
include mean transverse momentum fluctuations, correlation between transverse momentum 
and particle number, or even 2 and 3 particle correlation functions. This should be the subject 
of future work. 
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APPENDIX A: CONVERGENCE STUDY 

Not only for the sake of completeness we discuss in this section the convergence of various 
quantities with the sample size, i.e. the number of events, Ng.^ents, in out Monte Garlo scheme. 
Here we analyze final state (stable against electromagnetic and weak decays) particles only. We 
mainly take a closer look at the data sub-set of 20 ■ 2 ■ 10^ events, with A = Vi/Vg = 0.875 for 
the size of the bath, which already has been discussed in Section IVIIi 

There is a degree of freedom at so how to estimate the statistical uncertainty on the moments 
of a distribution of observables of a finite sample. The approach taken here is straight forward, 
but could, however, certainly be improved. 

In Fig. (fT9l) we show the evolution of the mean values (A^+) (left) and the variances ((AA^+)^) 
{right) of the distributions of positively charged hadrons for the 5 transverse momentum bins 
ApT,i, defined in Table [U with the sample size. Mean values of particle multiplicities in respec- 
tive bins are in rather good approximation equal to each other, but are, however, not identical 
due to finite resolution on the underlying momentum spectrum, even for A = 0.875 (bins were 
constructed using GGE events from an independent run). Variances converge steadily and 
are different in different bins, see Section IVIII The event output was iteratively stored in 
histograms, which were then evaluated after steps of 2 ■ 10^ events. 
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FIG. 19: Step histogram showing the convergence of the mean values (A^+) (left) and variances 
((AA^+)^) (right) for positively charged final state hadrons in transverse momentum bins Apx.i for a 
hadron resonance gas with A = Vi/Vg = 0.875. 
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FIG. 20: Step histogram showing the convergence of the scaled variance u;+ (left) of positively charged 

hadrons and the correlation coefficient between positively and negatively charged hadrons (right) 

in transverse momentum bins ApT.i for a final state hadron resonance gas with A = Vi/Vg = 0.875. 



In Fig. fl20|) we show the evolution of the scaled variance u;+ of positively charged final state 
particles (left) and the correlation coefficient p+_ between positively and negatively charged 
particles (right). The results for the respective transverse momentum bins can be compared to 
the second to last markers Figs. (11511161) . left panels, which denote the corresponding results of 
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FIG. 21: Histogram showing the results for the scaled variance a;+ {left) of positively charged hadrons 

and the correlation coefficient between positively and negatively charged hadrons (right) in the 

transverse momentum bin ApT^5 for a final state hadron resonance gas with A = Vi/Vg = 0.875. 200 
Monte Carlo runs of 2 • 10^ events each are analyzed. 

grouping the same data into 20 Monte Carlo sets of 2 • 10^ events each. 

In Fig. lpTl) we show the distribution of scaled variances of positively charged particles u;+ 
[left) and correlation coefficients between positively and negatively charged particles p_|__ (right), 
resulting from grouping again the same data set into 200 samples of 2 ■ 10^ events each. We 
chose the transverse momentum bin ApT,5 for a final state hadron resonance gas with A = 
Vi/Vg = 0.875. 

Monte Carlo results for A = 0.875 of the analysis shown in Fig. (l2T!) . are for the scaled variance 
u}+(ApT,5) = 0.8069 ± 0.0514, and the correlation coefficient p+_(ApT,5) = -0.0026 ± 0.0421. 
They are nicely scattered around the mean values, denoted by the bottom lines in Fig. fl20p . 
uj+(ApT,5) = 0.8082, and p+_(Apr,5) = -0.0028 respectively. 

They are also compatible with the analysis shown in Figs. fll5lll6p . of Section IVIIt 
u!+(ApT,5) = 0.8081 ± 0.0149, and p+_(ApT,5) = -0.0022 ± 0.0125, at the same value of A. 
The comparatively large statistical error on the analysis in Fig. fl2T|) is due to the splitting up 
into many small sub-samples. The mean values of different analyses agree rather well. 

Lastly, we show in Fig.( l22l) the results of additional Monte Carlo runs for values of A closer 
to unity. This time we have performed 20 runs of 1 ■ 10^ primordial events for A = 0.925, 0.950, 
and 0.975. As discussed above, error bars diverge, but convergence seems to be rather good. 
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FIG. 22: Evolution of the primordial scaled variance a;+ of positively charged hadrons (left) and the 

primordial correlation coefficient between positively and negatively charged hadrons {right) with 

the Monte Carlo parameter A = Vi/Vg in different rapidity bins Ayj. The solid lines show an analytic 
extrapolation from GCE results (A = 0) to the MCE limit (A — > 1). The 4 leftmost markers and their 
error bars represent the results of 20 Monte Carlo runs of 2 • 10^ events. 3 additional values of A have 
been investigated with 20 Monte Carlo runs of 1 • 10^ events. The rightmost markers denote the results 
of the extrapolation. 

The additional data has not been used for the extrapolation, so it can serve as an un-biased 
cross-check. 



APPENDIX B: THE CANONICAL BOLTZMANN GAS 



An analytical and instructive example is the canonical classical relativistic particle anti- 



particle gas discussed in 
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We use this example to show that, although the procedure 
is formally independent of one's choice of Lagrange multipliers, it is most efficient for those 
defined by Maxwell's relations. We start off with Eqs.([T]), and then discuss, in turn, the first 
and second moments of the multiplicity distribution of particles, and the first four moments of 
the Monte Carlo weight factor. 

The canonical partition function Zn-^^{Vi, P,Qi) of a system with volume Vi, temperature 
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T = (3 ^, charge Qi, particle number A^^i, and anti-particle number Mi = A'^i — Qi, is given by: 

The single particle partition function is given by Eq. (|3Tl) . tp = ^ m? j3~^ K2 {m(3). The canon- 
ical partition function with arbitrary particle number, but still fixed charge Qi, is obtained by: 

ZiV,,P,Qi) = Yl ZM,iV,,P,Qi) = lQ,{2V^ij) . (B2) 

JVi=Qi 

Here Iq^ is a modified Bessel function. Temperature is the same in both subsystems; the bath 
and the observable part. The partition function of the bath is therefore: 

Z{y2.(3,Q2) = lQA2V2i') . (B3) 

Imposing the constraints V2 = Vg — Vi, and Q2 = Qg — Qi, similar to Eq.(jl]), we find [39| for 

the canonical partition function, Eq.(l5]), of the combined system: 

00 

Z{Vg,P,Qg) = Yl hA'^V,^) /q^_q,(2(V,-Vi)^) = iQ^'^Vg^) , (B4) 

Qi=— 00 

as required. The weight factor is then: 

Analogous to Eq.(I7j) we find for the joint particle multiplicity and charge distribution: 



W{V^,Q,■,Vg,Qg\P) = ^^^^ . (B5) 



P{Qi,N,) = W{Vi,Qi;Vg,Qg\P) Zn,{VuP,Qi) . (B6) 



1. Monte Carlo Weight 

We next introduce Eq.( fT2l) . the joint GCE distribution of charges and particle multiplicity: 

PgceiQuN,) = ^^^^-^Z;v,(V^i,/5,gi). (B7) 
The Monte Carlo weight, Eq. fllSI) . is then given by: 

>VQi;Q.(y^.yj^^^) ^ W{V,,Q,;Vg,Qg\P) Z{V^,P,fi) e-'^^'^f . (B8) 
In accordance with Eq. lfTT]) . the distribution Eq.( 1B6l) is then equivalently written as: 

P{QuN,) = W'^'''^^{Vi;Vg\^i,(3) PgUQuNi) . (B9) 
The GCE partition function is: 

00 

Z(VS,A/x) = Yl e'^'^'' Z{Vi,P,Qi) = exp[yi2cosh(/3/i)] . (BIO) 
Qi=— 00 
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2. Moments of Distributions 



To define the multiplicity moments of the distributions Eq. flB6p or Eq. flBQP we write: 



oo oo 



w) = E E 

Ari=OQi=-oo 

Additionally we define the moments of the weight Eq. (1B5I) : 

oo oo 



(Bll) 



Ni=OQi=-oo 

and of the Monte Carlo weight Eq. (IBSP : 



(B12) 



oo oo 



(W") = J2 E [^'^''''^'{Vi;Vg\(3,f,)Y P,UQuN,) . (B13) 

Ni=OQi=-cxi 

We first attend to the first two moments of the multiplicity distribution. Substituting Eq. (lB6P 
or Eq.(lB9]) into Eq.(lBlI]) yields: 



and 



(N^,) = {V,^) '^-^ ^'^^^^ 



(B14) 



(B15) 



Canonical suppression of yields and fluctuations acts on the global volume Vg. In the GCE the 
first two moments are (Ni) = Viipe'^^, and (Nf) = (Ni)^ + (Ni), respectively. The CE limit is 
obtained by Vg Vi, and Qg = (Qi). Substituting Eq.([B14l) and Eq. flBTSj) into Eq.([52]), and 
using Eq.®, A = Vi/Vg, yields: 



UJ = \ UJce + (1 — A) Ul, 



gee 1 



(B16) 



where the CE scaled variance Uce of the combined system is given by 32|, 



UJr. 



1 



{2Vgij) 



{2Vgij) 



(BIT) 



and ujgce = 1 is the GCE scaled variance, as the particle number distribution is a Poissonian. 

We next apply our Monte Carlo scheme to an observable subsystem of volume Vi = 
embedded into a system of volume Vg = 75 /m^, charge Qg = 10, and temperature T = = 
O.lQOGeV . Particles and anti-particles have mass m = O.lAOGeV and degeneracy factor g = I. 
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The average charge content in the observable subsystem is then (Qi) ~ 6.667. The mean 
particle multiplicity, Eq.( lB14[) . is (A''i) ~ 7.335, and the scaled variance of particle number 
fluctuations, Eq. flB16l) . is u; ^ 0.3896. We will sample the GCE in Vi for various values of fiq 
and use the Monte Carlo weight, Eq. (]B8[) . to transform these samples to have the statistical 
properties required by Eq. flB6l) or Eq. flBQI) . For each value of /ig we have generated 50 samples 
of 2000 events each to allow for calculation of a statistical uncertainty estimate. 




FIG. 23: The first four moments of the Monte Carlo weight, Eq. ijBSp (left) and the first two moments 
of multiplicity distributions (right) , as described in the text. 



In Fig. fl23|) (right) we show, in open symbols, the mean value (A^i) and the variance ((AA^^i)^) 
of the particle multiplicity distribution of the original GCE samples for different values of 
chemical potential fig. The closed symbols denote mean value and variance of these samples 
after the transformation Eq. flBSP was applied. Independent of the original sample the result 
stays (within error bars) the same. However the statistical error is lowest for a chemical potential 
close to: 

i.e. when the initial sample is already similar (at least in terms of mean values) to the desired 
sample. This is reflected in the moments of the Monte Carlo Weight factor. Fig. (125]) (left). 
Higher moments have a strong minimum around fiq = 0.1896(76^, i.e. the weights are most 
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homogeneously distributed amongst events, and most efficient used is made of them. 
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